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We review the assumptions on which the Monte Carlo 
renormalization technique is based, in particular the analytic- 
ity of the block spin transformations. On this basis, we select 
an optimized Kadanoff blocking rule in combination with the 
simulation of a d = 3 Ising model with reduced corrections to 
scaling. This is achieved by including interactions with second 
and third neighbors. As a consequence of the improved ana- 
lyticity properties, this Monte Carlo renormalization method 
yields a fast convergence and a high accuracy. The results for 
the critical exponents are yn = 2.481(1) and yr = 1.585(3). 
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Applications of the Monte Carlo renormalization group 
(MCRG) to the three-dimensional Ising model Jl]-|j have 
become increasingly elaborate and complicated, and tend 
to require considerable computer resources. Neverthe- 
less, there are still uncertainties due to the basic as- 
sumptions underlying the renormalization transforma- 
tions used. In particular we will focus on the question 
concerning the analyticity of the transformation, which 
is closely related to the question whether the corrections 
to scaling vanish at the fixed point of the transformation. 
H In this letter we first calculate the analytic part of a 
divergent observable; this demonstrates that gross non- 
analyticities are normally absent. However, in general 
we may expect weak nonanalyticities due to corrections 
to scaling. Thus, second, we minimize their effects by 
adjusting the transformation as well as the Hamiltonian 
that is simulated. 

The MCRG method has amply been reviewed, fjljTj 
and here we only briefly outline the method. The reduced 
Hamiltonian is written as 



H(K , K 1 ,K 2r --;S) = -J2 K <*S a 



(1) 



where S is a spin configuration, the K a are couplings, 
and the S a are the conjugate lattice sums over spin prod- 
ucts, e.g. K\ is the magnetic field and Si = Sj the 
sum over all spins; K 2 is the nearest-neighbor coupling 
and S 2 = J2<nn> SiS j t ne sum over all nearest-neighbor 
pairs (si,Sj). A special 'coupling' is the background en- 
ergy density Kq; Sq is the number of spins. Application 
of a block-spin transformation to Monte Carlo generated 
configurations S leads to configurations S' described by a 
Hamiltonian W = H(K' Q , K[,K' 2 , • • • ; S'). The renormal- 
ized couplings K' a are assumed to be analytic functions 



of the original ones, even at the infinite system critical 
point. However, this property remains unproven in gen- 
eral, even for K' , or the so-called 'analytic part' of the 
transformation! 

It is straightforward to calculate, using Monte Carlo, 



B% = ((S' a S' )) = ((S' a - (S' a ))(S'e - {S' p )))) (2) 



and 



C% = ((S' a Sp)) = ((S' a (S' a )){S fi ~ (Sp))) (3) 



These lattice sum correlations are related || to the lin- 
earized transformation 



T Q/3 = dK'JdKp 



(4) 



(5) 



The dummy index summation rule applies to Greek in- 
dices. The matrix T is approximated by solving Eq. (|5|) 
after truncation to a finite number of couplings. Under it- 
eration of the block-spin transformation, the K a (a > 0) 
are assumed to approach a critical fixed-point, where the 
eigenvalues of T determine the critical exponents. 

Thus, the MCRG method relies on assumptions of 1) 
analyticity of the transformation; 2) convergence with 
the dimensionality of the coupling subspace; and 3) con- 
vergence to a critical fixed point. 

Concerning the third assumption, numerical work in- 
volving several subsequent transformations suggests 
that convergence to a fixed point does occur, and is de- 
scribed by an irrelevant exponent yi in the range —0.8 to 
-1.0. 

In order to investigate the second assumption, the 
number of couplings n c used in the analyses has in- 
creased considerably over the years; from 7 in Ref. 
to 99 in Ref. pj. A criterion to distinguish 'important' 
and 'less important' couplings was introduced in Ref. [||. 
The 'importance index' of an n-spin coupling is given by 
(2™/ 2 r where r is the average distance between the 
spins. This formula accounts for the facts that couplings 
tend to become less important when more spins are in- 
volved and when r increases. || An ordering according to 
this index leads to fast convergence || with increasing n c . 
On this basis we have restricted the present calculations 
to 20 even and 15 odd couplings, and indeed we observe 
good convergence with n c for all results presented here. 
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Next, we search for nonanalyticities in the 'analytic' 
part of the transformation when applied to an infinite 
system. First we investigate if the analytic parts of 
the susceptibility and the specific heat of the nearest- 
neighbor model, which are proportional to Ton and T022 
respectively are bounded at criticality. We express these 
quantities in derivatives of In Z and apply the chain rule: 

«9 2 lnZ d dlnZ d\nZ d 2 \nZ 

dK a 8Kp ~ 8K~ a 7/3 ~ T ^ a0 ~dK^ + T ^8K a dK^ 

(6) 

where T 7Q/ 3 = dT 1 ^jdK a . The derivatives of InZ can 
trivially be expressed in connected lattice sum correla- 
tions: 

{(S a S fj )) = T 7Q/3 (S;) + T 10 ((S a S' y )) (7) 

The T ia p are the only unknowns in Eq. (0); the corre- 
lations follow from the simulation, and the T 7j g from the 
standard MCRG analysis. There are not enough equa- 
tions to solve for To mm , but we can calculate the quantity 
A mm = L^ 3 T imm {S! y ) , in which the effect of a possible 
divergence of To mrra vanishes only in the case of unlikely 
cancellations. The factor L -3 normalizes S!y with respect 
to the system size L. 

We have done such calculations using the DISP, [^0|Jll| 
a special-purpose computer for Metropolis simulations 
of Ising models. The transformation is defined by 
the probability P(s') of a block spin s': P(s') — 
expos' Sb)/2cosh(ojSf,) where s& is the sum of the spins 
in a 2 3 block. It approaches the majority rule for large u>. 
In the limit of small u>, the block spins become indepen- 
dent and the critical singularity moves to the 'analytic' 
part. ]12| Numerical results for An and A22 did not sug- 
gest divergences in the analytic parts of the susceptibility, 
except where expected: for small w. Fig. 1 shows the nu- 
merical results for An, as a function of co for L = 8, 16 
and 32. 

A more strict test uses explicit calculation of Xb mm . 
We apply the chain rule to the second differentiation in 
d 3 In Z/dKadKpdK^ and express the derivatives of \nZ 
again in connected lattice sum correlations: 

((SaSpSy)) = TsapdSiySg)) + Ts/3((S a S' 7 S' s }} (8) 

Choosing a — [3 — m one can solve the unknowns Ts mm 
(5 > 0) from the numerical data, and thus isolate the 
term with To mm in Eq. ^. No signs of divergences are 
seen in the analytic part of the specific heat, except for 
small u, as illustrated in Fig. 2. 

These results are gratifying but the transformation 
may still have a weaker kind of singularity at the infinite 
system critical point. This can indeed be expected if the 
usual corrections to scaling are present at the fixed point 
of the block-spin transformation. As noted by Fisher 
and Randeria, || in general corrections to scaling van- 
ish only at the fixed point of an analytic transformation. 



Do they vanish at the fixed point of a block-spin trans- 
formation? This seems doubtful, in particular when the 
transformation contains free parameters which move the 
fixed point over the critical surface. Since the irrelevant 
fields are absent at a fixed point, any corrections should 
be due to some other mechanism. Weak singularities as- 
sociated with corrections to scaling can enter into the 
renormalized Hamiltonian via a weak nonanalyticity of 
the transformation. 

In order to suppress this problem in MCRG, we pro- 
pose the following strategy, which is applicable in a more 
general context than the Ising model: 

1. The Hamiltonian used to generate the Monte Carlo 
configurations is chosen such that the corrections to 
scaling are small. 

2. The transformation is chosen such that the fixed 
point is close to the Hamiltonian mentioned. 

To this purpose we included, in addition to nearest- 
neighbor couplings K2 = K nn , second and third neighbor 
couplings and K 3n in the Monte Carlo simulation, 
and optimized the ratio between the couplings, and the 
block-spin parameter uj. 

First we used a Monte Carlo Hamiltonian with K2 n 
and K 3n /K nn — 0.4, for which the corrections to scal- 
ing are small. Jl3|,[l4| Then the convergence to the fixed 
point, as apparent from the MCRG results for the eigen- 
values of T a p, becomes optimal for lo = 0.4. This is 
close to a variational value found by Kadanoff et al. 
The difference between the Monte Carlo and the fixed- 
point Hamiltonian follows, in a linear approximation, 
HH from the difference between (S a ) and (S' a ) as deter- 
mined from separate simulations of systems with com- 
patible sizes. This calculation was done in the coupling 
subspace (K nn , K 2n , K 3n ). 

A second approximation of the the u) — 0.4 fixed point 
was found by using a Monte Carlo Hamiltonian close 
to the first approximation. The fixed point was thus 
estimated (K nn , K 2n , K 3n ) = (0.1109,0.03308,0.01402). 
A finite-size scaling analysis of Monte Carlo results |f|] 
was used to determine the critical point more accurately: 
(K nn ,K 2n ,K 3n ) = (0.1114448,0.0332520,0.0140925), 
with a relative accuracy of 2 x 10~ 5 . This analysis showed 
that the corrections to scaling in the Binder cumulant |l6| ] 
are about 6 times smaller than for the nearest neighbor 
Hamiltonian. 

The bulk of the MCRG calculations took place at the 
estimated critical point, using system sizes L = 32, 16 
and 8 and lengths of 10 s , 2 x 10 7 and 10 7 sweeps re- 
spectively. The sensitivity to a variation in K nn was 
estimated from additional runs at K nn = 0.1114336 and 
K nn = 0.1114560. 

Further details, including the ordering of the cou- 
plings according to the importance index are contained 
in Ref. ||. Table 1 lists the resulting estimates for the 
exponents yx and yn, as determined from the largest 
eigenvalues of the T a p matrix. Statistical errors were 
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found by dividing the runs in 50 subruns. Finite-size 
and renormalization (approach of the fixed point) effects 
were determined with the procedures described e.g. in 
Rcf. |2j. The convergence of yj* and yn versus 2~" is, 
after correction for the finite-size effect, shown in Figs. 3 
and 4. For comparison we include results from Ref. ||, 
which used the standard nearest-neighbor Hamiltonian 
and the majority rule. Extrapolation of the data for 
L = 32 yields our final estimates for the critical expo- 
nents: y H = 2.481(1) and y T = 1.585(3). These results 
provide a significant improvement over previous MCRG 
analyses, not only concerning the statistical errors, but 
also the consistency with other results for the Ising uni- 
versality class which are summarized e.g. in Ref. [ fli) . 
The precise agreement of the renormalization results with 
those obtained by finite-size scaling confirms the validity 
of hyperscaling. |l7|| 

Earlier attempts to accelerate the convergence used op- 
timized transformations. fl8|-p3|| The present work uses 
this idea combined with a Hamiltonian with suppressed 
corrections to scaling. This alleviates the problem with 
the assumption of analyticity, and leads to a much im- 
proved convergence to the fixed point. The error due 
to the uncertainty margin of the exponent describing the 
renormalization effect practically vanishes. Furthermore, 
the rapid convergence to the fixed point eliminates the 
necessity of time-consuming simulations of large system 
sizes. Further improvements of the MCRG method may 
be possible by the introduction of more adjustable pa- 
rameters in the block-spin transformation, so that its 
fixed point can be moved to a point with even smaller 
corrections to scaling. 

This work is part of the research program of the 
"Stichting voor Fundamenteel Onderzoek der Materie 
(FOM)", which is financially supported by the "Ned- 
crlandse Organisatie voor Wetenschappelijk Onderzoek 
(NWO)" . 
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TABLE I. Numerical results for the renormalization expo- 
nents yT and yn, obtained after n block-spin transformations 
of a system of size L. 



exponent 


n 


L = 32 


L = 16 


L = 8 


yr 


1 


1.5885 (3) 


1.5885 (6) 


1.5868 (8) 


y T 


2 


1.5852 (5) 


1.5829 (9) 




yr 


3 


1.5829 (9) 






VH 


1 


2.48492 ( 4) 


2.48500 ( 7) 


2.48521 (22) 


yH 


2 


2.48309 (11) 


2.48327 (25) 




VH 


3 


2.48219 (27) 
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FIG. 1. The quantity An denned in the text versus the 1.65 
block-spin parameter ui for finite sizes L = 8 (□), L = 16 (A) 
and L — 32 (o). The lines are guides to the eye. Signs of a 
divergence with L appear only for small u. 




0.0 0.2 0.4 <x FIG. 4. The convergence of the magnetic exponent yu with 
m increasing number of iterations n of the block-spin transfor- 
mation. Results are shown for the present MCRG calculations 
(•) and for those reported in Ref. (o). 
FIG. 2. The quantity T022, which is proportional to the 2.50 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ 



analytic part of the specific heat, versus the block-spin pa- 
rameter W for finite sizes L = 8 (□), L — 16 (A) and L = 32 
(o). The lines are guides to the eye. Signs of a divergence 
with L appear only for small u. 
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FIG. 3. The convergence of the temperature exponent j/t 
with increasing number of iterations n of the block-spin trans- 
formation. Results are shown for the present MCRG calcula- 
tions (•) and for those reported in Ref. H (o). 
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